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Abstract 

We calculate the massive Wilson coefficients for the heavy flavor contributions to the 
non-singlet charged current deep-inelastic scattering structure function {x,Q'^) + 

xFl'- {x,Q‘^) in the asymptotic region ^ mf to 3-loop order in Quantum Chromo¬ 
dynamics (QCD) at general values of the Mellin variable N and the momentum fraction 
X. Besides the heavy quark pair production also the single heavy flavor excitation s ^ c 
contributes. Numerical results are presented for the charm quark contributions and conse¬ 
quences on the Gross-Llewellyn Smith sum rule are discussed. 



1 Introduction 


The non-singlet strnctnre fnnction xF^{x, Q^) has been a hrst experimental sonrce to determine 
the valence qnark distribntions Uv{x,Q^) and dy{x,Q‘^) of the nncleon in deep-inelastic charged 
cnrrent nentrino-nncleon scattering Flavor non-singlet distributions allow for clear mea¬ 

surements of the strong coupling constant = as/(47r) = g‘^/{4:7ry |^, due to the absence of 
gluonic effects in the QCD evolution [^. 

Massless and massive QCD corrections have been calculated in Refs. (6}]^ to hrst order in 
the strong coupling constant^ and in Refs. 


11-15 to 0(a^)[^ The massive O(a^) corrections 
were calculated in the asymptotic representation 16 , which is valid at high scales Q^. One may 


perform an 0{as) comparison for the process of single heavy quark excitation. Here the approx¬ 
imation holds for 17 , where = —q^ denotes the virtuality of the 4-momentum 


transfer g, and m is the heavy quark mass. 

The present data on the neutrino structure function xF^{x,Q‘^) or similar measurements at 
HERA have not yet reached the level of precision of 1-2 %, as is the case for the structure 
function F 2 {x,Q‘^) in neutral current deep-inelastic scattering. However, at neutrino factories 


planned for the future 19 this situation will change and even 3-loop QCD corrections will be of 


importance. The experimental precision reached for charged current interactions at the ep collider 
HERA can also be further rehned at future high energy facilities probing deep-inelastic scattering, 
such as LHeC 20 and EIC 21 . It is expected that non-singlet data taken at these facilities will 
help to improve further the knowledge of the strong coupling constant as(M|) [^. The 3-loop 
massless Wilson coefficient for the structure function combination xF^^ {x,Q‘^) + xF^ 
has been calculated in Ref. 23 . In the present paper we calculate the massive 3-loop Wilson 


coefficient for this combination in the asymptotic region ^ w?. 

The charged current structure functions xF^^ may be measured both in neutrino- 

and charged lepton-nucleon scattering. In the case of single gauge boson exchange the differential 
scattering cross sections read 
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(1.1) 
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( 1 , 2 ) 


X (l + (l-!,)2)fr(x,Q2) 
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Xr(x,Q^) ± (1 - (1 - !/)'0 XF,»'^(X,Q^) 


Here x = /ys and y = q.P/l.P denote the Bjorken variables, I and P are the incoming lepton 
and nucleon 4-momenta, and s = {I + P)^. Gp is the Fermi constant and M\y the mass of 
the hF-boson. F^"^{xyQ"^) are the structure functions, where the +(—) signs refer to incoming 
neutrinos (antineutrinos) and charged antileptons (leptons), respectively. We will consider the 
combination of structure functions 


xF, 


W+-W- I 


(x, QQ = xF^^ix, QQ + xF,^-{x, QQ 

were corrected in [^, see also [^. 
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^The massive 1-loop co rrec tions given in 
^Some results given in 13 have been corrected in Ref. 14 . 


( 1 , 3 ) 
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in the following. It can be measured projecting onto the kinematic factor hi = 1 — (1 — yY of 
the differential cross sections at x, = const, as a difference. It is given by 


pw+-w- 


{x.QY 


\Vdu?{d -d) + - 5 ) + Vu{u - u) 


„iv+-w-,Ns , TW+-w-,m 


\Vdc\Yd-d) + \Vsc\Ys-s) 


„W+-iy-,NS 

^g,3 




(1.4) 


with one massless Wilson coefficient ^ and two massive Wilson coefficients ^ 


and see Section 2, The coefficients Idj are the Cabibbo-Kobayashi-Maskawa (CKM) 


26,27 matrix elements, where 14 = |I4„p + and the present numerical values are [TS 


\Vdu\ 

= 0.97425 

(1.6) 

IKJ 

= 0.2253 

(1.6) 

IVdcl 

= 0.225 

(1.7) 

iKcl 

= 0.986 . 

(1.8) 


The massless Wilson coefficient reads to O(a^) 


C, 


VK+-W-,NS 

q,3 


"v 




(fc),lV+-W-,NS 

3 


k=l 


"v 


(1.9) 


where p denotes both the factorization and renormalization scale, which have been set equal 

/i = Pr = 


A{x) ® B{x) = / dxi / dx 26 {x — xiX 2 )A{xi)B{x 2 ) 


( 1 . 10 ) 


is the Mellin convolution, u, d and s are the light quark number densities and q denotes the 
corresponding antiquark densities, written here for a proton target. The valence distributions 
are given by 


qv = q-q ■ 

Very often one considers the case 

Si] 0 . 


(1.11) 

( 1 . 12 ) 


Various experiments are carried out using isoscalar targets, possessing neutron-proton symmetry. 
Here the following replacements have to be made 


u ^ \{u + d), u ^ \ {u + 

d —^ ^ (u -|- d), d —^ ^ (u dj^ . (1.13) 


The paper is organized as follows. In Section]^ we calculate the heavy flavor contributions to 
the non-singlet Wilson coefficients in the asymptotic region 3> m? to the structure function 
q 2 ^ 3 _ 1 qqp order in the strong coupling constant. We present the results both 
in Mellin N and x-space. Numerical results are given in Section Consequences for the 
Gross-Llewellyn sum rule are discussed in Section and Section contains the conclusions. An 
appendix deals with a technical detail. 


^For the scale dependence of the Wilson coefficient see e.g. 


28 
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2 The Wilson Coefficients 


The asymptotic heavy flavor corrections to the structure function ^ (x, are resulting 

from the Wilson coefficients and In the former case the charged current 

couples to a line of light quarks only, while in the latter case a heavy quark is excited by the 
exchanged W-boson with or without other heavy flavor effects through QCD corrections. In the 
following we work in Mellin space. The corresponding representation is obtained applying the 
transformation 

M[/(a:)](A^) = f dxx^~^f{x) (2.1) 

Jo 

to the above relations. One obtains 


+ 1 ) 
+ 1 ) 


4Ns ^m+-m-,NS 

^qq,Q^q,3 

4Ns ^m+-m-,NS 
^qq,Q'^q,3 


(N^ + 1) _ 


{Nf + 1) 


(Nf) 


( 2 . 2 ) 

(2.3) 


as the general relations. Here = A^^q{Nf + 1) denotes the massive non-singlet operator 

matrix element 29,30 . The expansion of both the massive operator matrix element (OME) and 
the massless Wilson coefficients in the strong coupling constant to 3-loop order leads to 






4(3).NS 

^qq,Q 


(2),NS^(l),W+-W-,NS/ , 1 N , A(3),m+-m-,NS 


^^qq,Q ^q,3 


1 -I- ttsCg ^. 


{nf +1)+ 




-fa. 


-fa„ 


^qq,Q + ^q,3 F + -*-) 


q2),m+-m-,Ns^ 
0,3 


.(3),NS .(2),NS^(l),W-^-W-,NS.,. ,x 
^qq,Q +^qq,Q ^ q,3 OVi^ + l) 


TC*, 


(3),m+-m-,NS. 


= L 


q,3 i^F + 1) 

m+-m-,NS 
q,3 


{Nf + 1) + Og 3 


W+-W-,NS 


(Nf), 


■ym+-m-,NS 
0,3 
convention 


. (2.4) 


( 2 . 6 ) 

( 2 . 6 ) 


with the massless Wilson coefficient up to 3-loop order. Here we use the 


/{Nf) = /{Nf + 1) - /{Nf) 


(2.7) 


The notation label ‘iVi? -f 1’ in some of the above quantities means that they depend on Nf 
massless and one massive flavor. 

The calculation of the different contributions to the Wilson coefficients (2.4, 2.5) is per¬ 
formed m D = A + e dimensions to regulate the Feynman integrals. The treatment of 75 has 
to be considered. In the flavor non-singlet case both for the massive OMEs and the massless 
Wilson coefficients, 75 always appears in traces along one quark line. For the asymptotic Wilson 
coefficient this line is massless, while for also heavy quarks contribute due to single 
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charged current flavor excitations like s —>■ c. Still we are considering the case in which power 
corrections k ^ N, k > 1, are disregarded and therefore the corresponding line has to 

be treated as massless. In this case a (reversed) Ward-Takahashi identity implies anti-commuting 
75 , mapping the corresponding vertex correction to a self-energy correction, which is the same 
in the case without 75 . 

The inclusive massive OME to 3-loop order for even and odd moments N has been cal- 
The corresponding Feynman integrals have been reduced using integration- 
31 applying an extension of the package Reduze 2 3^ The master integrals 


culated in Ref. 30 


by-parts relations 


have been calculated using hypergeometric, Mellin-Barnes and differential equation techniques, 
mapping them to recurrences, which have been solved by modern summation technologies using 
extensively the packages Sigma 35,36, EvaluateMultiSums, SumProduction 37, psum (^, 


and HarmonicSums 39 


In Mellin N space the Wilson coefficient can be expressed by nested harmonic sums Ss{N) 
which are defined by 


N 


Sb,a{N) = = 1, 6, Ui e Z, 6,Ui ^ 0, iV > 0, iV e N. 


k=l 


kW 


( 2 , 8 ) 


In the following, we drop the argument N of the harmonic sums and use the short-hand notation 
Ss{N) = Ss- The Wilson coefficients depend on the logarithms 


ro = >« (|r 


and Lm = In ( ^ 


(2.9) 


As a short-hand notation we define the leading order splitting function 7 gg^ without its color 
factor 

3 A ^2 + 3jv + 2 ' 


7, 


(0) ^ 4 


m 


2 Si - 


2N{N + 1) 


The massive Wilson coefficient for the structure function xK 


W+-W- 


( 2 . 10 ) 


(x, Q^) in the asymptotic 


region in the on-shell scheme in Mellin space is given by 

*(/V) = 

4(29iV2 + 29N - 6 ) 
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+ aU ClTp 


-*-^,( 0 ) 2 ,- (j 2 
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,128 64^ 
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^The package Reduze 2 uses the packages Fermat 
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and Ginac 
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where C 3’^ i^fp) denotes the contribution due to the massless Wilson coefficient at 3 -loop 


order 23 , and the polynomials P* are 


Pi 

= -17iV^-34Ar3-29Ar2-12iV-24 

( 2 . 12 ) 

P 2 

= -3Ar4 - 6iV3 - 47iV2 - 20iV + 12 

(2.13) 

7^3 

= Ar4 + 2 iV 3 - Ar 2 _ 2 Ar - 4 

(2.14) 

Pa 

= 3iV^ + 6Ar3 + 47iV2 + 20iV - 12 

(2.15) 

P 3 

= 7iV^ + 14Ar3 + 3iV2 - 4Ar - 4 

(2.16) 

Pe 

= 19A^^ + 38N^ - 9N‘^ - 20N + 4 

(2.17) 

P 7 

= 28Ar^ + 56iV^ + 28iV2 + 2 Ar + 1 

(2.18) 

Ps 

= 33iV^ + 54iV^ + 9iV2 - 52iV - 28 

(2.19) 

P 9 

= 57A^^ + 96iV^ + 65iV2 - lOiV - 24 

( 2 . 20 ) 

Pio 

= 112Ar^ + 224iV^ + 121iV2 + 9iV + 9 

( 2 . 21 ) 

Pii 

= 141Ar4 + 246iV3 + 241iV2 - SA^ - 84 

( 2 . 22 ) 

P 12 

= 181A^^ + 266iV^ + 82N^ - 3iV + 18 

(2.23) 

Pl3 

= 235Ar^ + 524iV^ + 211iV2 + 30iV + 72 

(2.24) 

Pl4 

= 359Ar^ + 772N^ + 335N^ + 30iV + 72 

(2.25) 

Pl5 

= 501A^^ + 894iV^ + 541iV2 - 116iV - 204 

(2.26) 
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Pi 6 = 561iV^ + 1122iV3 + 767iV2 + 302iV + 48 (2.27) 

Pi 7 = 1131iV^ + 2118iV3 + 1307iV2 + 32iV- 276 (2.28) 

Pi 8 = 1139iV^ + 2710iV3 + 635iV2 + 216iV + 828 (2.29) 

Pi 9 = 1199iV^ + 2398iV3 + 1181iV2 + 18iV + 90 (2.30) 

P 20 = 1220iV^ + 2359iV^ + 1934iV2 + 357iV - 138 (2.31) 

P 21 = 3iV^ + lliV^ + 10iV3 + 19iV2 + 23iV + 16 (2.32) 

P 22 = 12iV® + 16iV^ + 18iV3 - ISiV^ - 5iV - 8 (2.33) 

P 23 = 27iV® + 863iV^ + 1573iV^ + 1151iV2 + 144iV- 36 (2.34) 

P 24 = 648iV® - 2103iV^ - 4278iV=^ - 3505iV2 - 682iV - 432 (2.35) 

P 25 = -11145iV® - 32355iV^ - 37523iV^ - 14329iV^ + 1240iV2 - 1032iV - 2088 (2.36) 

P 26 = -151iV® - 469iV® - 181iV^ + 305iV3 + 80iV2 - 88N - 56 (2.37) 

P 27 = 7V6 + 3iv5 - 8 iV^ - 21iV3 - 23iV2 - 12iV - 4 (2.38) 

P 28 = 6iV® + 18iV® - iV^ - 20iV3 + 46Ar2 + 29iV - 6 (2.39) 

P 29 = 15Ar® + 36iV® + 30iV^ - 24iV3 + 3Ar2 + 16iV + 20 ^^2.40) 

P 30 = 155iV® + 465iV^ + 465iV^ + 155iV^ + 108iV2 + 108iV + 54 (2.41) 

P 31 = 216iV*^ + 567iV^ + 687iV^ + 285iV3 + 37iV2 - 44iV + 12 (2.42) 

P 32 = 309iV® + 807iV^ + 693iV^ - 463iV3 - 638iV2 + 68N + 216 (2.43) 

P 33 = 525iV® + 1575iV® + 1535iV^ + 973iV^ + 536iV2 + 48iV - 72 (2.44) 

P 34 = 609iV‘^ + 1485iV® + 1393iV^ + 83iV3 - 422iV2 + 156iV + 216 (2.45) 

P 35 = 795 iv 6 + 2043iV® + 2075iV^ + 517iV3 - 298N^ + 156iV + 216 (2.46) 

P 36 = 868 iV® + 2469iV® + 2487iV^ + 940iV^ + 171iV2 + 207iV + 144 (2.47) 

P 37 = 1407iV® + 3825iV® + 4211iV^ + 1783iV^ - 250iV2 - 240iV + 144 (2.48) 

P 38 = 1770iV® + 4671iV® + 4765iV^ + 1205iV3 - 227iV2 + 1044iV + 756 (2.49) 

P 39 = 7531iV® + 23673iV^ + 23055iV^ + 7375iV^ + 1614iV2 + 936iV - 324 (2.50) 

P 40 = -4785iV® - 19140iV^ - 18754iV® + 1320iV® + 12723iV^ + 6548iV3 + 4080iV2 

-648iV - 1728 (2.51) 

P 41 = -45iV® - 162iV^ - 858iV® - 936iV® - 1629iV^ - 1094iV3 - 804iV2 

-40iV + 192 (2.52) 

P 42 = 7V8 + 4iVV 13iV® + 25iV® + 57iV^ + 77iV=^ + 55iV2 + 20iV + 4 (2.53) 

P 43 = 3iV® + 12iVV 16iV® + 6 iV^ + 30iV^ + 64iV3 + 73iV2 + 40iV + 12 (2.54) 

P 44 = 3549iV® + 14196iV^ + 23870iV® + 25380iV® + 15165iV^ + 1712iV3 - 2016iV2 

+ 144iV + 432 (2.55) 

P 45 = 5487iV® + 21948iV^ + 36370iV® + 28836iV® + 11943iV^ + 4312iV^ + 2016iV2 

-144iV-432 (2.56) 

P 46 = 10807iV® + 43228iV^ + 63222iV® + 40150iV® + 14587iV^ + 9018iV^ + 7452iV2 

+2376iV + 324 (2.57) 

P 47 = 42591iV® + 166764iV^ + 245664iV® + 129982iV® - 13295iV^ - 25978iV3 

+3560iV2 - 3192iV - 4464 (2.58) 

P 48 = -18351iVi° - 89784Ar® - 210021iV® - 271638iV^ - 219369iV® - 90572iV^ 

-26491iV^ - 7790iV3 - 1992iV2 - 2760iV - 2160 (2.59) 
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F 49 = + 825iV^ + 109664iV® + 331682iV^ + 457641Ar® + 346145iV® 

+219290iV^ + 86724iV^ + 13608iV^ + 14256iV + 10368 (2.60) 

F 50 = S28N^° + 3492iV^ + 4305iV® - 2013iV^ - 8540A^® - 3S22N^ - 1157N^ 

-3057iV^ - 4:112N‘^ - 324:N + 576 (2.61) 

F 51 = 8274iV^° + 37149A^® + 53630A^® + 7538A^^ - 59902iV® - 55159iV'^ 

-6994iV^ + 3272N^ - 904SN^ - 1656A^ + 2160 . (2.62) 

Here Ca = Cp = {N^ - 1)/{2 Nc),Tf = 1/2 and = {N^ - l)(Ar/ - 4:)/N^ denote 

the color factors for SU{Nc), with A'"c = 3 in case of QCD. The constant F 4 is given by 


2 13 

F 4 = - 4(2 ln^(2) + - ln^(2) - —(4 + 1614 ( 7 


(2.63) 


which is related to mnltiple zeta valnes |^, 7 = l/^”) n G N, n > 2, denote the valnes 

of Riemann’s (/-fnnction at integer argnment and Li„(a:), n G N, is the polylogarithm 42 . Terms 
containing the color factor /N^ stem only from the massless Wilson coefficient and the 

non-singlet valence anomalous dimension at 3-loop order 43 , cf. Appendix A. 


44 45 . These expansions can be obtained 


One obtains the analytic continuation of the harmonic sums to complex values of N by 
performing their asymptotic expansion analytically, cf. 
automatically using the package HarmonicSums 
obey the shift relations 


39 . Furthermore, the nested harmonic sums 


SUN) = SUN - 1 ) + 


(sign(&)) ^ 

NW 


-S,{N) 


(2.64) 


through which any regular point in the complex plane can be reached using the analytic asymp¬ 
totic representation as input. The poles of the nested harmonic sums Ss{N) are located at the 
non-positive integers. In data analyses, one may thus encode the QCD evolution together with 
the Wilson coefficient for complex values of N analytically and hnally perform one numerical 
contour integral around the singularities of the problem]^ 

over 


In x-space the Wilson coefficient is represented in terms of harmonic polylogarithms 47 


the alphabet {/o, /i, /-iji which were again reduced applying the shuffle relations 48 
defined by 


They are 


Hb,a{x) = J dyfb{y)Hs{y), Ho^{x) = —U{x), = I, 

k 

fo{x) = fl{x) = f-iU = 

X 1 — X 


1 -1- X 


(2.65) 

( 2 . 66 ) 


The Wilson coefficient is represented by three contributions, the (...)+-distribution, the 5(1 — x)- 
term, and the regular term, with the -|—distribution being defined by 


It reads 



dy[F{y)]^g{y) 


dyF{y) [g{y) - g{l)] 


(2.67) 


^For precise numerical implementations of the analytic continuation of harmonic sums see 


46 
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( 2 . 68 ) 


Again, we used the short hand notation H^{x) = also here. The transformation of the Wilson 
coefficient to the MS scheme for the heavy quark mass affects the massive OME at 3-loops and 


was given in Ref. 30 ; the terms are the same in the unpolarized and polarized case. 


The -|—distribution of the Wilson coefficient also contains 1/(1 — x)2-terms, cf. 30 . The 


explicit expressions for the Wilson coefficient H^§{Nf + 1) can be easily obtained from Eq. (2.6) 


3 Numerical Results 


In what follows, we will choose the factorization and renormalization scale = Q^. We hrst 
study the behaviour of the massive and massless Wilson coefficients in the small and large x 
region and then give numerical illustrations in the whole x-region. 

At small X, the pure massive Wilson coefficient behaves like 


LfsiNp + 1 ) 


ans 

'^q,3 


(Np) oc ag-C^Tp ln 2 (a;)- 1 -a; 
o 




ln"^(x). 


(3.1) 


while in the region x —)■ 1 one obtain^ 

Lf,{Np + l)-C^^I{Np) oc a^CpTp 

+a^gCpTF 


224 80 82 ' 

IT + 9'^" + 3 " 


448 160 , 


1 — X 

ln 2 (l — x) 
1 — X 


®There is a typo in the second contribution to the 0{al) term of Eq. (3.2) in 49 , which should be oc ln((3^/m^). 


19 


























(3.2) 


Likewise one obtains for + 1) at small x 


Hf^liNF + l)-C^^^UNF + l) oc ai^CFTF\n\x) + ^ 


■yNS 


16 


^CaCfTf - fC^Tf \n\x), 


(3.3) 


and in the region a; —)■ 1 

H^I{Nf + 1)-C^^I{Nf + 1) oc alCFTF 

+a^gC‘^TF 


224 80 , 8 , 2 ' 

- H“ ^ A/T 

27 9 " ^ 3 ^ 


448 160 , _,2 


1 — X 

ln^(l — x) 
1 — X 


(3.4) 


The above results can be compared with the case of the massless Wilson coefficient (7^ 3 and 


CoS, respectively. The following dominant terms are found in the small and large x regions, 


oc al^CFTF\n\x) 

O JabcJabc 

CIs’^^\Nf) oc — —ln^(a;) 


2 rn i„2/ 


T5 Nr 


cII’^^\Nf) oc -a,2C'F ln(a;) 


<x a'i 


^NS,( 3 )/ 7 y ] ^ rr^ 
^q,3 K^^F) OC 


Icl - 2CaCf 


29 


ln^(a;) 


53 


'abc Jabc 


~r ~^p — 


15 


30 


2 

15 Nr 


-N, 


ln^(a;) 


and 


oc ai-CFTF 


oc aF-^Ci-Tr 


9 


ln^(l — x) 

1 — X 
ln"‘(l — x) 
1 — X 


oc a.4CF 

(iVr) « aFd 

C'HNi.NF) « aFCl 


ln(l — x) 

1 — X 
ln^(l — x) 
1 — X 
ln®(l — x) 
1 — X 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 


respectively. 

Let us hnally also consider the case Lq 7 ^ 0. Terms of this order have the following leading 
small and large x behaviour 


C^J'^‘^\Nf) oc -ai^-^CFTFLQ\n{x) 


16 


(3.15) 
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,.3 (^Vf) 
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'ln^(l — x) 


oc a'l^QClLq 


X 


respectively, and are less singular than those for Lq = 0. 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 



Figure 1: The corrections up to 3-loop order to the combination of the structure functions {x, Q^) -|- 

xFr {x, Q"^) off a proton target, including both the massless and the charm contributions in the asymptotic 
approximation in the on-shell scheme for rric = 1.59 GeV as a function of x and Q^. The parton 
distribution functions of Ref. 52 have been used with as{M^ = 0.1132. These settings are the same in 
the subsequent Figures. 


The small x behaviour can be compared with leading order predictions for the non-singlet 
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Figure 2: Ratio of the structure functions {x, + xF^ {x, off a proton target up to 3-loop 

order for the charm contribution and the massless terms for three flavors. 



Figure 3: Ratio of the structure functions xF^^{x,Q‘^) + xF^ {x,Q‘^) off a proton target up to 3- 
loop order for the charm contribution and the massless terms for three flavors from tree-level to 0(af) at 
= 100 GeV2. 


evolution kernel in Refs. 50,51 


principle pattern 


Cfco; 


fc-i-i 


In 


. Indeed, both the massive and massless contributions follow the 
(x). However, as is well known 51 , less singular terms widely 
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cancel the numerical effect of these leading terms. For the large x terms the massless terms 
exhibit a stronger soft singularity than the massive ones. 



Figure 4: The ratio of the structure functions {x, ix^) + xF^ to xF^^ {x,Q‘^) + 

xF^ (x,Q‘^), varying G [Q^/4,4Q^] at 0(af) and = 100 GeV^ as a function of x. 



Figure 5: The corrections up to 3-loop order to the combination of the structure functions xF^^ (x, Q^) + 
xF^ (x,Q^) off a nucleon in an isoscalar target, including both the massless and charm contributions in 
the asymptotic approximation. 
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In Figure 1 we present the combination of structure functions xF, 


W+-W-,NS 


up to O(a^) at a 


proton target as a function of x and up to 3-loop order for three massless flavors and charm 
in the asymptotic representation, assuming rric = 1.59 GeV 53 in the on-shell scheme for mass 
renormalization and setting fF = Q^. It shows valence-like scaling violations moving towards 
smaller values of x. 

The effect of the charm contributions are illustrated in comparison to the purely massless 
ones up to 0{al) in Figure 2 evolving with Q^. The contribution amounts to a correction of up 
to +3% in the small x region and to —3% in the large x region with some evolution in Q^. 

In Figure 3 we show a similar ratio as in Figure 2, but taken for each order of the perturbative 
corrections to xF^^ {x, -|- xF^ separately at = 100 GeV^. In the small x region 

the corrections grow from 1% to 3%. At large values of x the relative corrections tend to become 
negative and do also amount to a few per cent. Figure 4 shows the remaining scale variations of 
the combination xF^"^ [x, + xF^ up to 3-loop order in the region /i^ G [Q^/4,4Q^] 

normalized to the value at /i^ = (yellow band) as a function of x for = 100 GeV^. 
The behaviour is overall flat with ±1% variations in the medium x-range and shows larger 
uncertainties at very low and large values of x. The behaviour is similar at other virtualities. 

Figure 5 shows the combination of structure functions xF^^~^ up to O(a^) at a nucleon 
in an isoscalar target. The differences to the case of the proton target shown in Figure 1 turn 
out to be rather small in the case of this combination, which widely pronounces isoscalarity due 
to the present combination of currents, if Sy ~ 0. 


4 The Gross-Llewellyn Smith Sum Rule 


The Gross-Llewellyn Smith sum rule 54 refers to the hrst moment of the flavor non-singlet 
combination 


dx [F 3 "^(x, Q2 ) + F^/{x, g2)] = 6GGLs(h.), 


(4.1) 


with ds = Us/ti and idealized GKM mixing. The 1-loop [^[^|^[^, 2-loop 57 , 3-loop 58 and 
4-loop QGD corrections 1^ ,60 in the massless case are given by 


GGLs(a.), = l-ds + a2(-4.58333 + 0.33333Nf) + a^(-41.4399 + 7.74370A^i7 - 0.17747A^^) 
+a^(-479.448 + 140.796A^f - 8.39702Ar2 ^ o.l0374Ar3) , (4.2) 

choosing the renormalization scale /i^ = Q^, for SU{3)c. Here Np denotes the number of active 
light flavors. The expression for general color factors was given in Ref. 59,60. Note that 


the QGD corrections to the Gross-Llewellyn Smith sum rule and to the polarized Bjorken sum 
rule 61 are identical up to 0(af). 


For the asymptotic massive corrections (2.4, 2.5) only the hrst moment of the massless Wilson 


coefficients contributes, weighted by the hrst moments of the valence quark densities 

and the corresponding GKM matrix elements, since the hrst moments of the massive non-singlet 
OMEs vanish due to fermion number conservation. This also holds at higher order. Unlike the 
case of the polarized Bjorken sum rule 
contribute. One obtains 


49 , in the present case two massive Wilson coefficients 


dx [F^^ix,Q‘^) + F^^ix,Q^)] = 2[(|Urf,|2 + \Vsu?){uy) + (|U,/ + \Vdu?){dy) 

xGGLs(hs, Np + 1) 


(4.3) 
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= 2 [2 ■ 0.9999 + 0.9998] CGLs(as, Np + 1) 


(4.4) 


for the charm excitation with additional Np light flavors. Here we assumed = 0 and 
(g„) = Jg dxq^{x) denotes the hrst moment of a valence distribution. The CKM matrix ele¬ 
ments currently imply a very small deviation from the factor 6 in (4.4) in the asymptotic case 
for proton targets. Very similar results are obtained for isoscalar targets. There are, in 
particular, no logarithmic corrections in the mass scales involved. The heavy quark just induces 
a shift Np —)• Np + 1 w.r.t. the case of Np light quarks only. Different results are obtained in 


the non-inclusive tagged-flavor case 13 , which has been calculated to O(a^). Here no inclusive 
structure functions are considered unlike the case in Ref. 
in the tagged-flavor case were derived in Refs. 


62,63 


14 . Corresponding power corrections 


5 Conclusions 

We calculated the heavy flavor non-singlet Wilson coefficients for the charged current structure 
function xF^^~^ {x, Q^) to O(a^) in the asymptotic region m^. They can be expressed in 

terms of nested harmonic sums and poly logarithms. In contrast to the neutral current case, here, 
a second heavy flavor Wilson coefficient contributes, which describes the flavor excitation due 
to a s —)■ c transition. The heavy flavor contributions to this combination of structure functions 
amounts up to the 0(±3%) for proton targets. Very similar results are obtained for deuteron 
targets, as the present combination approximately selects isoscalarity. In the small and large x 
regions the heavy flavor effects at 3-loop order are visible on top of those up to 2-loop orders 
by an enhancement and a depletion of 0(1%), respectively, in the ratio to the 3-loop massless 
Wilson coefficient. This is below the present resolution reached at HERA, but may become of 
importance in high luminosity measurements at future colliders. The scale dependence reached 
for the non-singlet combination studied turns out to be widely flat. We also presented the leading 
terms in the heavy and light flavor Wilson coefficients at small and large values of the momentum 
fraction x. 

In the asymptotic case the contribution of the massive Wilson coefficients to 

the Gross-Llewellyn Smith sum rule reduce to the replacement Np —)■ Np -|- 1 in the massless 
approximation. One has to note CKM-matrix effects here. 


A Appendix 

The massless Wilson coefficients and anomalous dimensions for xF^^ {x,Q‘^) -|- xF^^ {x,Q‘^) 
contain a new color factor proportional to (= 40/3 in QCD) [^|^. This color factor 

does not appear in the massive operator matrix elements at 3-loop order, which we would like 
to discuss in the following. Those terms, however, do arise in individual diagrams contributing 
to the OMEs but they cancel in the complete result. The cancellation can be seen as follows. 

The color factor can appear when there are two separate fermion lines connected by three 
gluons. One fermion line connected to three gluons yields the color structure 

Tr[nV] = . (A.l) 

Thus, two fermion lines connected by three individual gluon propagators produce the color 
structure 

^abc jabc 



Tr [FfF] ^Tr [F fF] = 
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For each such diagram there is a corresponding diagram with the fermion flow along the closed 
fermion loop reversed. An example of a pair of diagrams is given in Figure The closed loop 
has to have three quark-gluon vertices in order to produce the color factor and therefore the loop 
must have three fermion propagators. Keeping the direction of the momenta hxed, the reversal 
of the fermion flow entails a change of the sign of the momentum in the numerator of each 
fermion propagator, 


pj — m? pj — m? 


(A.3) 



Figure 6: Example of a pair of diagrams which each contain a term proportional to the color 
factor 


This sign can be factored out since traces over an odd number of Dirac matrices vanish and 
yields a global factor (—1). Besides that, the reversal of the fermion flow also reverses the order 
of the color generators A in Eq. (A.2) which flips the sign in front of the term, but 


leaves the term unchanged. We see that each pair of such diagrams has exactly the same 

integrand, but the sign in front of the color factor is changed. Therefore, this color factor 

cancels in the sum for each pair of diagrams. 
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